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Abstract 


This paper reports on research for extending the Cooperative Con- 
trol Synthesis methodology to include a more accurate modeling of the 
pilot’s controller dynamics. Cooperative Control Synthesis (CCS) 
is a methodology that addresses the problem of how to design con- 
trol laws for piloted , high-order , multivariate systems and/or non- 
conventional dynamic configurations in the absence of flying qualities 
specifications. This is accomplished by emphasizing the parallel struc- 
ture inherent in any pilot- controlled, augmented vehicle . The original 
CCS methodology is extended to include the modified optimal control 
model (MOCM), which is based upon the optimal control model of 
the human operator developed by Kleinman, Baron , and Levison in 
1970. This model provides a modeling of the pilot’s compensation 
dynamics that is more accurate than the simplified pilot dynamic 
representation currently in the CCS methodology. Inclusion of the 
MOCM into the CCS also enables the modeling of pilot- observation 
perception thresholds and pilot- observation attention allocation ef- 
fects. This Extended Cooperative Control Synthesis (ECCS) allows 
for the direct calculation of pilot and system open - and closed-loop 
transfer functions in pole/zero form and is readily implemented in 
current software capable of analysis and design for dynamic systems. 
Example results based upon synthesizing an augmentation control law 
for an acceleration command system in a compensatory tracking task 
using the ECCS are compared with a similar synthesis performed by 
using the original CCS methodology. The ECCS is shown to provide 
augmentation control laws that yield more favorable, predicted closed- 
loop flying qualities and tracking performance than those synthesized 
using the original CCS methodology. 


Introduction 

Increasing aircraft agility, maneuvering at high angles of attack, and exploring radical flight 
vehicle geometries to obtain low observability are areas of current research that show promise of 
greatly increasing aircraft mission performance. To fully exploit these and other possible new 
capabilities, future aircraft may require high-order, multivariate flight control systems or demand 
designs dealing with nonconventional flight dynamics. Although design guidance is available in 
the form of flying qualities specifications for aircraft exhibiting conventional dynamics, very little 
design guidance is available to the flight control designer for synthesizing control laws to achieve 
both good piloted performance and good flying qualities for high-order and/or nonconventional 
configurations. 

In 1979, Schmidt proposed a synthesis methodology that addresses the problem of how to 
design control laws for piloted, high-order, multivariate and/or nonconventional dynamics con- 
figurations in the absence of flying qualities specifications. This methodology, referred to as 
“Cooperative Control Synthesis (CCS),” emphasizes the parallel structure inherent in any pilot- 
controlled augmented vehicle. (See fig. 1.) The CCS methodology is applicable to high-order 
systems and leads to control laws for good piloted performance and subjective evaluation. In 
this method, optimal control theory is utilized for both the control law synthesis and a simplified 
modeling of the pilot’s compensation dynamics. By including this simplified model, the CCS 
methodology explicitly includes design objectives based upon pilot acceptability. The original 
CCS methodology was extended by Innocenti and Schmidt (1984) to include state estimation in 




Figure 1. Conceptual block diagram of Cooperative Control Synthesis. 


the pilot model and to allow a measurement feedback control law. Garg and Schmidt (1989) have 
shown how the CCS methodology can be used to synthesize pilot-acceptable display dynamics. 

This paper reports on research for extending the CCS methodology to include a more accurate 
modeling of the pilot’s controller dynamics than the simplified pilot dynamic representation 
currently used in the CCS. This is accomplished by replacing the simplified model of the pilot’s 
controller dynamics currently in the CCS with the modified optimal control model (MOCM) 
(Davidson and Schmidt, 1992). The MOCM is based upon the optimal control model (OCM) of 
the human operator developed by Kleinman, Baron, and Levison (1970) and Baron, Kleinman 
and Levison (1970). Inclusion of the MOCM into the CCS also enables the modeling of 
pilot-observation perception thresholds and pilot-observation attention allocation effects. This 
extended CCS allows for the direct calculation of pilot and system transfer functions in pole/zero 

form and is designed for easy implementation in current software capable of analysis and design 
for dynamic systems. 

A theoretical development of the Extended Cooperative Control Synthesis (ECCS) method- 
ology is provided, and this methodology is used to synthesize augmentation control laws for an 
acceleration command system in a compensatory tracking task. This analysis is compared with 
similar designs performed by using the original CCS methodology, and conclusions are presented. 

Symbols and Abbreviations 

system dynamic matrix 
system control matrix 
system output matrix 
system disturbance matrix 
tracking error 
Kalman filter gain matrix 
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cost function control-rate weighting 

gain matrix 

ith pilot regulator gain 

objective function 

regulator Ricatti solution matrix 

transfer function gain 

augmented pilot control gain vector 

pilot control gain vector 

system measurement matrix 

control vector dimension 

disturbance vector dimension 

state vector dimension 

output vector dimension 

measurement vector dimension 

cost function state weighting matrix 

augmented weighting matrix 

augmented state vector 

cost function control weighting 

cost function augmentation control weighting 

Laplace variable 

time 

pilot output 

pilot-commanded control 
motor noise intensity matrix 
observation noise intensity matrix 
pilot disturbance vector 
motor noise disturbance 
observation noise disturbance vector 
state noise intensity matrix 
augmented noise intensity matrix 
state disturbance vector 
plant and disturbance state vector 
Pade delay state vector 
pilot observation vector 
measurement vector 
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6 input to plant 

6 pitch attitude 

p signal- to- noise ratio 

5] filter Ricatti solution matrix 

r effective time delay 

Tf) “neuromotor” lag (neuro-lag) 

X augmented state vector 

Abbreviations: 

CCS Cooperative Control Synthesis 

Eqq expected value 

ECCS Extended Cooperative Control Synthesis 

LQG linear quadratic Gaussian 

MOCM modified optimal control model 

OCM optimal control model 

PR pilot rating 

rms root mean square 


Subscripts: 

a augmented 

c pilot-commanded control 

cmd command signal 

d delay 

obs pilot observed 

p pilot 

5 plant and delay augmented system 

u control 

S plant input 

0 control-rate augmented system 

1 plant and pilot augmented system 

Operators and superscripts: 

T transpose 

— 1 inverse 

* optimal 

(.) derivative with respect to time 

(.) estimate 
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Theoretical Development 

The Cooperative Control Synthesis (CCS) methodology involves the simultaneous solution 
of two coupled optimal control problems (Papavassilopoulos, Medanic, and Cruz, 1979). One 
optimal controller can be thought of as representing a pilot’s control dynamics, whereas the 
other represents the augmentation control law dynamics. (See fig. 1.) A simultaneous solution is 
required because the pilot’s control strategy is a function of the augmented vehicle dynamics, and 
the vehicle augmentation control law iff not known a priori. This section presents a development 
of the Extended Cooperative Control Synthesis (ECCS) methodology incorporating the modified 
optimal control model (MOCM) to represent both the pilot’s control dynamics and a direct- 
output-feedback linear quadratic controller for the augmentation control law dynamics. 


Control Solution 


The plant dynamics to be acted upon by the two optimal controllers acting in parallel, 
augmented with the system disturbance dynamics, are given by the state space, time-invariant 
linear form 


x = Ax + B<5 + Ew ’ 

6 = 6p + 6 a 

y = Cx > 

yobs = y + v y 

z = Mx , 


( 1 ) 


where x is an n^-dimensional state vector, 6 is an n^-dimensional vector equal to the sum of the 
pilot’s control input ( 6 p ) and augmentation controller input (6 a ), and w is an n w -dimensional 
disturbance vector modeled as zero-mean Gaussian white noise with an intensity W. The vector 
y of dimension n y represents variables that the pilot can perceive, either by observation or feel. 
The perceptual model observed by the pilot (y Q bs) assumed to be corrupted by an observation 
noise (v y ), i.e., a zero-mean Gaussian white noise process with intensity V y . The vector z of 
dimension n z denotes system measurements available for feedback. For this development, the 
assumption is made that the system is stabilizable and detectable and that the measurements 
are noise free. The basic structural components of these two controllers acting on the dynamic 
system are represented in the block diagram in figure 1. 

Solution for pilot dynamics. This formulation of the CCS methodology incorporates the 
MOCM (Davidson and Schmidt, 1992) to model the pilot’s compensation dynamics. A block 
diagram of the model components of the MOCM is given in figure 2. The MOCM is based upon 
the optimal control model (OCM) of the human operator developed by Kleinman, Baron, and 
Levison (1970) and is a variation of simplified optimal pilot models developed by Hess (1976), 
Broussard and Stengel (1977), and Schmidt (1979). This model of the human operator is input 
compatible with the OCM and retains other key aspects of the OCM, such as the modeling of 
pilot-observation perception thresholds and pilot-observation attention allocation effects. Unlike 
the OCM, however, the structure allows for the direct calculation of pilot and system transfer 
functions in pole-zero form. 

In the MOCM, the pilot’s effective time delay is modeled by a second-order Pade approxima- 
tion. This time delay is placed at each of the pilot’s outputs and is treated as part of the plant 
dynamics for the determination of the pilot’s regulation and filter gains. To simplify the notation, 
this development considers the case of a single control input, although the algorithm can easily 
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Figure 2. Pilot’s block diagram in modified optimal control model. 


be extended to account for multiple pilot and/or augmentation control inputs. A second-order 
Pade approximation is given by 


UJ _ 1 - 1/2 (ra) + 1/8 (ts) 2 
Up 1 + 1/2 (ts) + l/8(rs) 2 

where r is the delay interval, u p is the pilot’s output, and u d is the delayed pilot’s output. In 
state space form, this can be expressed as 

x d = A d x d + B d u p 
6 P = u d = C d x. d + up 

where is a two-element vector of Pade delay states. 

The plant dynamics augmented with the pilot’s effective time delay are given by 
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z = [M 0]{ x *} 


The pilot’s observation vector is given by 


(4) 


Yobs 




+ 


This model makes the assumption that the pilot’s control task can be defined by the 
minimization of the quadratic performance index ( J p ) given by 


J P = Eoo {y r Q P y + UpTpUp + iipfpUp | 


(5) 
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subject to pilot observations (y G bs) with cost-functional weightings Qp > 0 , r p > 0, and f p > 0. 
The system given by equation (4) can be expressed in a control-rate formulation as 


d_ 

dt 
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Yobs = [C 0 0] { x d } + 


ILr 


z = [M 0 0] < 
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or by defining a new state vector as 


X = 


T T 
x x d u p 


leads to 

X = AqX -j- B pii p + B a <5 a + E 0 W ^ 
yobs = Cl* + Vy > (6) 

z = MjX J 

The pilot controls the augmented aircraft. The augmentation control law is given by 


6 a = G a z = G a MiX 


(7) 


The system with the augmentation loop closed is given by 

X = (Aq + B a G a Mi )X + B p ii p + E 0 w > 

Yobs = ClX + Vy 

z = MiX * 


(8) 


The performance index can be rewritten in terms of the augmented state vector (X) as 


where 


Jp — Eqo ^ X QX + Up fpiip ^ 


Q = diag 


Cf QpCi, r p 


(9) 


The minimizing control law is obtained by application of standard LQG solution techniques 
(Kwakernaak and Sivan, 1972) to the delayed augmented system (eq. (8)). This leads to the 
linear feedback law (Kleinman, Baron, and Levison, 1970) 


n 

Up — — 9n+\ u p — ^ ^ 

i= 1 


( 10 ) 
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where n = n x + 2 (system states plus two Pade states) and X is the estimate of the state X. By 
letting 


1 


Tn = 


and 

one obtains 


9n + 1 

Ip — r f)[9\ y • • • > <7n] 


T V Up + Up = U c 

where the pilot’s commanded control ( u c ) is given by 


U c ipX s 


( 11 ) 

( 12 ) 


C T TV 

i x dJ ' 


where x s • 

The n + 1 feedback gains (G p ) are obtained from 

J fjl 

= [$1 j • • ■ » 9ni 9n- |-l] = fp Bp 
where K is the unique positive semidefinite solution of the Ricatti equation 


0 = (A 0 + B a G a M,) J K + K(A 0 + B a G a Mi) + Q- KBp/^BjK 


(13) 


(14) 


To account for the uncertainty of the human operator’s control input, motor noise (v u ) is added 
to the commanded control (u c ). Thus, 

Tyjilp T Up ~ Uc + Vu (1^) 

where v u is a zero-mean Gaussian white noise process with intensity Y u . Solving for ii p gives 

(16) 


-1 1 1 

u p = — Up H u c + — v u 




•V 


The controller gains are assumed not to be affected by the inclusion of the motor noise. This 
reduces the solution of the human operator model to a subopt imal control law. 


Therefore, the augmented system combined with equation (16) is given by 
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or 

X = (Ai + B a G a Mj)X + B \u c + EjWi 
The current estimate of the state (X) is given by a Kalman filter 

X = (Ai + B a G a Mi)X 4- Bi u c 4- F(y 0 bs — y) 

X = (A^ 4- B a G a Mi — FCi)X 4* FC^X 4* B \u c 4~ Fv^ 

F = ECfV " 1 


(17) 


(18) 


where 
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The covariance matrix of the estimation error (E) is the unique positive semidefinite solution 
of the Ricatti equation 

0 = (Ai + B a G a M!)E + E(Ai + B a G a M,) T + EjWjEf - ECfV^CiE (19) 
where Wj = diag[W, V u ] with W > 0, V u > 0, and V y > 0. 

Solution for augmentation controller. As in the original Cooperative Control Synthesis 
methodology, this formulation employs a direct-output-feedback linear quadratic controller 
(Levine and Athans, 1970) for the augmentation controller. This augmentation control law 
is chosen to be optimal with respect to the objective function ( J a ) given by 

Ja = Roc {y Qay d" Up T a \Up + 6 a r a 2 ^a + faUp^ (20) 

where Q a > 0, r a \ > 0, r a 2 > 0, and f a > 0. This approach hypothesizes that the pilot’s index 
of performance is correlated with the pilot’s subjective rating of the vehicle handling qualities 
in the design task (Schmidt, 1981). This can be accomplished by letting Q a = Q p , r al = r p , 
and f a = f p \ then, J a can be expressed by 


Ja — Jp d- Eqq ■|<5jV a 2^a^ (21) 

The solution for the augmentation controller is carried out on the system in the presence of the 
pilot’s control compensation. The system with the augmentation loop open is given by 
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X — AjX d B a £ a d Biu c + E 1 W 1 (22) 

The pilot’s compensation dynamics are given by 

X = (Aj + B„G a Mi)X + Bju c + F(y 0 b s — y) 
u c = ~l\X 


where l \ = [l p 0]. Therefore, the pilot-augmented system is given by 


d_ 

dt 



Ai -B1/1 

FCi Ai + BaGaMi-Bjii-FCj 
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or, by letting 


then 


q = [* T x T ] T 

q — Aq + B<5 a + Ew 
z = Mq 


(23) 
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The augmentation control objective function written in terms of the augmented state vector (q) 
is given as 

(24) 


Ja = Eoo jq r Qq + 6^r a2 6 a } 


where 


Q = 


CiQpCf 0 0 

0 r v 0 


0 


0 G^f a G 


The minimizing control law is given by 

6 a = G q z = G n Mq 

where the minimizing gains (G a ) are given by 

G„ = -r a2 1 B r HLM r (MLM T )' 1 


(25) 


(26) 


The matrices H and L are obtained from the simultaneous solution of two coupled steady-state 
Lyapunov equations: 


(A + BG a M) L + L (A + BG a M) T + E W E T = 0 
(A + BG a M) T H + H (A + BG a M) + (q + M T G^r a G a M) = 0 


(27) 


System closed-loop and pilot state space models . A state space representation of the 
closed-loop pilot-vehicle system is given by 




Aj + B a G a Mi — BjZi 

FCi Ai + B q G q M] — BiL — FCj 


ira-r: :h::i 


Cj 0 

Cg OJ 




(28) 


where 


Cg = [G a M C d 1 0 0 0] 


A block diagram of the model components of the pilot’s dynamics is given in figure 2, and a 
state space representation of the pilot’s dynamics is given by 
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( 29 ) 
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where 



Algorithm Implementation 

The solution to the ECCS problem requires the simultaneous solution for the pilot’s control 
dynamics (by solution of the MOCM) and the augmentation controller (by solving a direct- 
output-feedback control synthesis problem). The solution of the MOCM involves solving two 
algebraic Ricatti equations for the pilot’s control and estimation gains. Solution of the direct- 
out put- feedback, linear quadratic control synthesis problem requires the simultaneous solution 
of two Lyapunov equations. The optimal augmentation gains can be obtained by numerical 
iteration by employing a conjugate gradient search technique (Fletcher and Powell, 1963). A 
conceptual flowchart of the ECCS algorithm is given in the appendix. The implementation of 
this algorithm requires the capability of solving steady-state Ricatti and Lyapunov equations. 

The structure of the extended CCS methodology allows implementation in current software 
capable of analysis and design for dynamic systems. Implementation in this type of environment 
allows for rapid calculation of pilot and system transfer function descriptions from state space 
models and rapid determination of system frequency responses. Also, this environment allows 
users to interactively modify various pilot and plant parameters and quickly ascertain the impact 
of these changes on the pilot/closed-loop performance. 

Values of effective time delay, “neuromotor” lag (neuro-lag), observation, and motor noise 
intensities are chosen in the same manner as for the OCM. Like in the OCM, desired values 
of neuromotor lag (r T? ) can be obtained by an appropriate choice of cost-function control-rate 
weighting ( f p ). Manual control experiments have shown that the effective time delay of the pilot 
(r) is typically 0.1 to 0.2 sec (Kleinman, Baron, and Levison, 1970). 

The intensity of the observation noise (' V y ) is dependent upon the nature of the display, 
human limitations, and the pilot’s environment. Over a wide range of viewing conditions, the 
diagonal elements of the observation noise-intensity matrix are proportional to the variance of 
its associated observed variable. Single- axis manual- tracking control tasks have shown that on 
the average, p Vi = 0.01, which corresponds to a normalized observation noise ratio of —20 dB. 
The intensity of the motor noise (V u ) is assumed to be proportional to the variance of the 
commanded control (u c ). An analysis of single-axis, manual-tracking control task experiments 
has shown that, typically, p u = 0.003, which corresponds to a normalized motor noise ratio of 
—25 dB (Kleinman, Baron, and Levison, 1970). 

The next section will present an application of the ECCS method to synthesize pilot optimal 
control gains for an acceleration command system and a comparison of these results with results 
obtained by using the original CCS methodology. 

Synthesis Example 

This section will present the application of the ECCS methodology to synthesize “pilot 
optimal” augmentation control laws for an acceleration command system in a compensatory 
tracking task. The ECCS results will also be compared with a similar synthesis performed by 
using the original CCS methodology. 


Example 


For an acceleration command system, the controlled system dynamics in transfer function 
form is given by 


e 

6 


k_ 


(30) 
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To be consistent with the references (Schmidt, 1979, and Innocenti and Schmidt, 1984) the 
plant gain is set to k = 11.7. ’ 

In this task, the pilot’s objective is to track a displayed command signal. The signal to be 
tracked is generated by a second-order, low-pass filter driven by unit intensity white noise 

^cmd + 3.O0 cmd + 2.250 cmd = 3.67 w(t) (31) 

The state space representation of the plant augmented with the disturbance dynamics is 
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where 

^ j^crndi ^cmdi 

and 6 is the sum of the pilot’s and the augmentation controller inputs, 8 = 8p + 8 a . 

The pilot’s performance index is given by 

Jp — Eoo |(# cm d — 0) 2 + fp8p | (33) 

and the augmentation control performance index is given by 

Ja = Jp + Eqo |r a 2^ 2 | (34) 

In this development, the common assumption is made that the human controller is able to 
perceive both position error and error rate from the display. Therefore, the pilot’s observation 
vector is given by 


_ / ^cmd ^ \ f e \ [10— 1 01 

^ Un,d-Sj ”“U/ +V, ' = [o 1 0 -lJ X + 


(35) 


The measurements available for feedback are z = {0,()} T . Therefore, the augmentation 
control law is 


6 a - G a z = Gq | q | = G a Mx = G a 


0 0 10 
0 0 0 1 


(36) 


where 


— [Gfl G^] 

Values of effective time delay, neuromuscular time constant, and observation and motor noise 
intensities were chosen to be consistent with single- axis manual-tracking control tasks. The 
control-rate weighting on 6 p in J p , which is f p , was adjusted to obtain a neuromuscular time 
constant of t v = 0.1 sec. The pilot’s effective time delay (r) was chosen to be 0.1 sec, and the 
pi ot s observation noise-to-signal ratio and motor noise-to-signal ratio were chosen to be —20 dB 
and -25 dB, respectively. Values of input parameters for the MOCM are summarized in table 1. 

Results obtained from applying the ECCS methodology to this system are given in table 2 
Synthesis results from applying the original CCS methodology to this system are given in table 3. 
Model-based predictions of pilot root-mean-square (rms) performance, augmentation gains, and 
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augmented plant poles are shown as a function of augmentation control weighting. The predicted 
rms tracking error and pilot control are determined from the steady-state covariance of the closed- 
loop pilot-vehicle system. Also shown is a prediction of the subjective pilot rating (PR), based 
upon the value of the pilot’s index of performance (Schmidt, 1979), calculated using the relation 

PR = 25 ln(10 Jp) + 0.3 (37) 

These predicted subjective pilot ratings are based upon values of the pilot’s index of performance 
calculated using both the ECCS methodology in table 1 and the original CCS methodology in 
table 2. 


Table 1. MOCM Input Parameters 


Effective time delay, (r), sec 0.1 

Neuromotor lag, (r n ), sec 0.1 

Observation noise ratio, dB —20 

Motor noise ratio, dB —25 

System disturbance (0 C mt \/w) 3.67/(s 2 + 3s 4- 2.25) 

Disturbance intensity (W) 1.0 

Objective function observation weights (Q p ) diag[l, 0] 

Objective function input weight (r p ) 0 


Table 2. ECCS Augmentation Results for Acceleration Command System 


Control 

weighting 

rms 

tracking 

error 

rms 

pilot 

control 

j p 

Predicted 

PR 

G e 

G e 

Augmented 

poles 

(a) 

Unaugmented 

0.8222 

0.8474 

0.9363 

5.9 

0 

0 

(0),(0) 

100.0 

.7807 

.7991 

.8368 

5.6 

-.0073 

-.0043 

[0.29, 0.09] 

10.0 

.6104 

.6181 

.4894 

4.3 

-.0464 

-.0346 

[0.74,0.27] 

1.0 

.4938 

.5636 

.3095 

3.1 

-.1148 

-.1179 

[1.16,0.59] 

.1 

.4463 

.7089 

.2529 

2.6 

-.2890 

-.2680 

[1.76,0.92] 

.01 

.4301 

.8736 

.2379 

2.4 

-.4102 

-.3877 

(2.20), (2.86) 


a (-) = Stable real pole; [•,•] = [u\£] = FYequency and damping, respectively, of complex pole pair. 


Table 3. Original CCS Augmentation Results for Acceleration Command System 
[Data from Innocenti and Schmidt (1984)] 


Control 

weighting 

rms 

tracking 

error 

rms 

pilot 

control 

Jp 

Predicted 

PR 

Go 

G 0 

Augmented 

poles 

(a) 

Unaugmented 

1.17 

1.0 

1.86 

8 

0 

0 

(0),(0) 

6 100.0 








10.0 

.79 

.66 

.77 

5 

-.06 

-.03 

[0.9,0.19] 

1.0 

.63 

.49 

.46 

4 

-.22 

-.08 

[1.6,0.30] 

.1 

.62 

.61 

.45 

4 

-.58 

-.13 

[2.6,0.30] 

.01 

.61 

.80 

.44 

4 

-.79 

-.15 

[3.0,0.30] 


a ( ) = Stable real pole; [•, ■] = [w, £] = FYequency and damping, respectively, of complex pole pair. 
^No data available. 
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A plot of predicted pilot rating versus augmentation gain magnitude for various values of 
augmentation control weighting (r a 2 ) is given in figure 3. Figure 4 shows a locus of closed-loop 
system eigenvalues as a function of augmentation control weighting. Frequency responses of 
the unaugmented and augmented plant, the pilot compensation dynamics, and the system loop 
transfer functions as a function of augmentation control weighting (r a 2 ) are given in figures 5-7. 

Discussion of Example 

An analysis of the results obtained with CCS incorporating the MOCM (table 2) shows an 
improvement both in terms of predicted rms tracking performance and in predicted subjective 
rating over those obtained from the original CCS formulation (table 3). For example, a synthesis 
design that achieves a predicted pilot rating of Cooper-Harper level I (Cooper and Harper, 1969) 
at a minimum of gain magnitude is at a control weighting of 1.0. (See fig. 3.) At this design 
point, the predicted average Cooper-Harper pilot rating is a 3 and the rms tracking performance 
is 0.4938, which compares with a Cooper-Harper rating of 4 and an rms tracking performance 
of 0.63 with the original CCS methodology. 
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Figure 3. Predicted pilot rating versus gain 
magnitude. 


Figure 4. Augmented system eigenvalues as a 
function of control weighting. 
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Figure 5. Augmented system frequency response as a function of control weighting. 
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Figure 6. Pilot compensation frequency response as a function of control weighting. 


r al 

— 0.01 
— 1 
— 100 



Figure 7. System loop frequency response as a function of control weighting. 


As can be seen from the augmented system frequency responses (fig. 5), the augmented plant 
dynamics are type 0. The pilot compensation required to achieve the k/s loop shape in the 
region of crossover (approximately 3 rad/sec) is given in figure 6. 

Research in pilot-in-the-loop systems by McRuer (1980) has shown that in a compensatory 
tracking task, a human will adjust his dynamic compensation such that the system loop dynamics 
will approach k/s in the region of frequency crossover. When this is achieved, the system has 
good closed-loop stability characteristics. The system loop frequency responses (fig. 7) show that 
this has been achieved in all the designs for the values of control augmentation weight chosen. 

Concluding Remarks 

This paper has presented research for extending the Cooperative Control Synthesis methodol- 
ogy to include a more accurate modeling of the pilot’s controller dynamics. Cooperative Control 
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Synthesis (CCS) is a methodology that addresses the problem of how to design control laws 
for piloted, high-order, multivariate systems and/or nonconventional dynamic configurations in 
the absence of flying qualities specifications by emphasizing the parallel structure inherent in 
any pilot-controlled augmented vehicle. The simplified model of the pilot’s controller dynamics 
currently in the CCS is replaced by the modified optimal control model (MOCM). The MOCM 
(based upon the optimal control model of the human operator developed by Kleinman, Baron, 
and Levison in 1970) provides a modeling of the pilot’s compensation dynamics that is more 
accurate than the simplified pilot dynamic representation currently in the CCS methodology. 
Inclusion of the MOCM into the CCS also enables the modeling of pilot-observation perception 
thresholds and pilot observation attention allocation effects. 

The structure of this Extended Cooperative Control Synthesis (ECCS) methodology allows 
implementation in current software capable of analysis and design for dynamic systems. Im- 
plementation in this type of environment allows for the rapid calculation of pilot and system 
transfer function descriptions from state space models and the rapid determination of system 
frequency responses. Also, this environment allows users to modify interactively various pilot 
and plant parameters and to quickly ascertain the impact of these changes on the closed-loop 
pilot-vehicle performance. 

The ECCS methodology was used to synthesize pilot optimal augmentation control laws for 
a simple dynamic example in a compensatory tracking task. This analysis is compared with 
similar designs using the original CCS methodology. Analysis results obtained with the CCS 
incorporating the MOCM show an improvement both in terms of predicted root-mean-square 
tracking performance and predicted subjective rating over those obtained from the original CCS 
formulation. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
March 4, 1994 
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Appendix 

Methodology Flowchart of Extended Cooperative Control Synthesis 

A conceptual flowchart of the ECCS algorithm is presented in this appendix. The imple- 
mentation of this algorithm requires the capability of solving steady-state Ricatti and Lyapunov 
equations. 
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